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Neutralino dark matter is well motivated, but also sufTers from two shortcomings: it requires 
gravity-mediated supersymmetry breaking, which generically violates flavor constraints, and its 
thermal relic density fl is typically too large. We propose a simple solution to both problems: neu- 
tralinos freezeout with Q ~ 10—100, but then decay to --^ 1 GeV gravitinos, which are simultaneously 
light enough to satisfy flavor constraints and heavy enough to be all of dark matter. This scenario is 
naturally realized in high-scale gauge-mediation models, ameliorates small scale structure problems, 
and implies that "cosmologically excluded" models may, in fact, be cosmologically preferred. 
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Supersymmetric extensions of the standard model of 
particle physics are among the prime candidates for new 
microphysics. Among their many virtues, supersymmet- 
ric models naturally predict new particles that are candi- 
dates for dark matter. The most well studied of these are 
thermal relic ncutralinos [l[ , superpartners of the Higgs 
and electroweak gauge bosons. The thermal relic density 
of ncutralinos is dependent on unknown supersymme- 
try parameters. However, order-of-magnitude estimates 
yield relic densities that are consistent with ^ 



f^DM/i = 0.1050 



+0.0041 
0.0040 



(la) 



(1) 



where JIdm is the observed energy density of non- 
baryonic dark matter in units of the critical density, and 
h ~ 0.73 is the normalized Hubble parameter. This re- 
markable fact has not only motivated supersymmetry, 
but has also focused attention on "cosmologically pre- 
ferred" models, in which the neutralino thermal relic den- 
sity is exactly that required for dark matter. Such studies 
have implications for a large range of experiments, from 
direct and indirect dark matter searches to those at the 
Large Hadron Colhder (LHC) at CERN. 

The neutralino dark matter scenario is not without 
its blemishes, however. First, for the neutraUno to be 
stable, it must be the lightest supersymmetric particle 
(LSP). In particular, it must be lighter than the grav- 
itino. This requires gravity-mediated supersymmetry 
breaking models, in which low energy bounds on fla- 
vor and CP violation are generically violated by several 
orders of magnitude. Gauge-mediated supersymmetry 
breaking (GMSB) models l3| elegantly avoid these con- 
straints, but such models have gravitino LSPs and so are 
incompatible with neutralino dark matter. 

Second, although general arguments imply that the 
neutralino thermal relic density is of the right order of 
magnitude, in concrete models, it is often too large: 
Ncutralinos are Majorana fermions, and so annihila- 
tion to quarks and leptons is P-wave suppressed. In 
addition, gauge coupling unification and radiative elec- 
troweak symmetry breaking typically imply Bino-like 



ncutralinos, which suppresses annihilation to gauge and 
Higgs bosons. These effects together enhance relic densi- 
ties to values that may far exceed those given in Eq. ([T]) . 

These two shortcomings of neutralino dark matter are 
usually considered unrelated and addressed separately. 
One may, for example, consider gravity-mediated sce- 
narios, such as minimal supergravity, where low energy 
constraints are satisfied by unification assumptions. One 
then further focuses on special regions of parameter space 
in which the neutralino relic density is reduced to ac- 
ceptable levels through, for example, resonant annihi- 
lation [4], stau co-annihilation or significant Bino- 
Higgsino mixing [9] . Alternatively, one may simply aban- 
don the hope that the order-of-magnitude correctness of 
the neutralino thermal relic density is a significant lead 
in the hunt for dark matter and explore other mecha- 
nisms for dark matter production. For example, one may 
consider GMSB models with thermally produced graviti- 
nos 7]. (Note, however, that recent Lyman-a constraints 
requiring > 2 keV Q imply that the gravitino ther- 
mal relic density O^/i^ « 1.2 (m^/keV) must be signifi- 
cantly diluted through late entropy production 9] for this 
possibility to be viable.) More recently, GMSB- like mod- 
els with gravitino dark matter produced by late decaying 
gauge singlets have also been proposed [10l|. 

In this work, we consider the possibility that the two 
shortcomings described above are not separate issues, but 
are in fact pointing to a single resolution. We propose 
that ncutralinos do, in fact, freezeout with very large 
densities. However, they then decay to gravitinos, which 
are light enough to accommodate the GMSB solution to 
the flavor and CP problems, but heavy enough to be all 
of dark matter. In analogy to Goldilocks planets, which 
have temperatures that lie within the narrow window re- 
quired to support life, these supersymmetric models have 
gravitino masses in the narrow window required to satisfy 
both particle physics and cosmological constrants, and so 
we call this "Goldilocks Supersymmetry." 

The essential features of this scenario may be illus- 
trated by simple scaling arguments. Consider models 



2 



in which there are two mass scales: the scale of the 
standard model superpartner masses to, and the grav- 
itino mass m^,. The freezeout density of neutralinos 
is inversely proportional to the neutralino annihilation 
cross section, and so by dimensional analysis, fl^h^ ~ 
(ctw)^^ ~ TO^. The gravitino relic density is therefore 
flijh^ = {mQ/'fh)^l^h? ^ m^rh. At the same time, a 
natural solution to the supersymmetric flavor and CP 
problems requires to^ <C rh. We find, then, that we can 
always make large enough to explain dark matter by 
raising to^, and to together with their ratio fixed. The 
essential question, then, is whether the scenario may be 
realized with to TeV, as required for a natural solution 
to the gauge hierarchy problem, and whether it passes all 
other particle physics and astrophysical constraints. 

To analyze this question concretely, we consider the 
example of minimal GMSB models ll| . Such models are 



specified by the 4+1 parameters M^, A, iVm, tan/3, and 
sign(^), where Mm is the messenger mass, A — F/M^ti, 
where F is the supersymmetry breaking scale in the mes- 
senger sector, iVm is the number of 5 + 5 messenger pairs, 
tan/3 = (77")/(_ff°), and fj. is the supersymmetric Hig- 
gsino mass. In terms of these parameters, the gauge- 
mediated contributions to squark and slepton masses are 



167r2 



where C( — f^^, with hypercharge Y = 
C- = for gauge singlets, | for SU(2) l doublets, and | 
for SU(3)c' triplets. The gaugino masses are 



(2) 



Tjj, and 



167r2 



(3) 



where i = 1,2,3 for the U(l)y, SU(2)l, and SU(3)c 
groups, ci = |, and C2 = cs — 1. As indicated, these 
masses are generated at the energy scale Mm. We deter- 
mine physical masses through renormalization group evo- 
lution to the weak scale and radiative electroweak sym- 
metry breaking with SoftSUSY 2.0 

In addition to the gauge-mediated masses, there are 
gravity-mediated contributions. These generate the grav- 

where Fq is the total super- 



Fq 

V3A/. 



itmo mass to^ 

symmetry breaking scale and ~ 2.4 x 10^* GeV is 
the reduced Planck mass. Because Fq receives contribu- 
tions from all supersymmetry breaking i^-terms, Fq > F. 
For direct gauge mediation, Fq ~ F, but this is model- 
dependent. Here, we assume Fq — F, and so 



F 



M„,A 



V3M, V3Af* 



(4) 



Our results are not changed significantly for _Fo > F. 

Gravity-mediation also generates flavor- and CP- 
violating squark and slepton mass parameters (to{^ )ab, 
where i,j = 1,2,3 label generation, A,B = L,R label 



chirality, and f = l,u, d. The chirality-preserving param- 
eters are naturally ~ to^; for concreteness, we assume 
\{m{j)LL\ = \{m{j)RR\ = m^. The chirality-violating 
masses require the breaking of electroweak gauge sym- 
metry (and possibly horizontal symmetries); we assume 
|(to^)l/j| ^ XfjiTiQ, where the A^- are Yukawa couplings. 
Finally, we assume 0(1) CP-violating phases for both the 
gravity- and gauge-mediated masses, as detailed below. 

Given these assumptions, the most stringent con- 
straints are the flavor-changing observables Am^ and 
ex, and the CP- violating, but flavor-preserving, electron 
and neutron electric dipole moments (EDMs) 13.,.14,.15|: 



Am^j}^^^ < 3.5 X 10"^^ MeV 

g|USY ^ 2.3 X 10"^ 

4 < 1.6 X 10"^^ e cm 

dn < 2.9 X 10"^^ e cm . 



(5) 
(6) 
(7) 
(8) 



In the mass insertion approximation, these constrain 
{SL)ab = {"miAAB/frif, where rhf is an average / 



mass 



The leading constraints are from AniK on 
Re[(<5f2)LL((5f2)flfl], from on lm[iSl^)LLiSf^)RR], 
and from the EDMs on the gauge-mediated masses. 

The supersymmetric contributions to the kaon ob- 
servables are Ato|!^^"^ = Re(M) and e|I^^^ = 
Im(M)/(v^ATO^P), with M as given in Ref. ;i6]. For 
concreteness, we choose the 6 phases to maximize the su- 
persymmetric contribution for each kaon observable. The 
constraints from Attik and ex are therefore not simul- 
taneously applicable, but the most stringent constraint 
smoothly interpolates between these as the phase varies. 
For the EDMs, we first use micrOMEGAs 1.3.7 [l3| to 



determine the supersymmetric contribution to a^, the 
anomalous magnetic moment of the muon. The EDMs 
are, then, de = j^a^teLTi9cp and d„ = |(4(irf + (i„), 
where dd and du are determined from dg with a — > as. 
Ml — > M3, TOj ^ TOi ~, and the introduction of appropri- 
ate color factors [la]. We set tan^cp = 1 hi the EDMs. 
Note that the EDMs may be suppressed, depending, for 
example, on the origin of the ^ and B parameters. 

The resulting constraints are given in Fig. [T] The ob- 
servables AniK and ex require to^ < 30 GeV (500 GeV) 
for neutralino mass ~ 100 GeV (1 TeV). In contrast, 
the EDMs are insensitive to TOq, since they do not rely 
on gravity-mediated contributions. They are found to re- 
quire to^ ^ 1 TeV, in agreement with earlier work [l^ . 
These results are, of course, subject to the assumptions 
we have made. However, they imply that in any model 
in which gravity-mediated contributions are at their nat- 
ural scale and all mass parameters have 0(1) phases, the 
standard model superpartners must be heavy, and the 
LSP is the gravitino, not the neutralino. 

For Am = 1, the lightest standard model superpartner 
is the lightest neutralino x throughout parameter space. 
In Fig. [1] we also show the freezeout density flyh^, that 
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FIG. 1: Neutralino thermal relic density 0.^1? in the (m^, A) 
plane, for iVm = 1, tan/3 = 10, > and mt = 175 GeV. The 
right-hand axis gives the neutralino mass ~ 1.3 x 10~^A. 
Regions to the right of the ea' and Am/f contours and below 
the de and dn contours are disfavored. The neutralino is the 
LSP in the shaded region. 



is, the relic density if neutralinos were stable, determined 
using micrOMEGAs 17| . These results illustrate the dif- 
ficulties for neutralino dark matter. At the weak scale, 
typically /i, > Mi, and x is Bino-like. Its annihilation 
is therefore suppressed for the reasons noted above. For 
m-)^ = 100 GeV, 17^ /i^ ~ 1 is already far too large, and 



for the heavier superpartner masses favored by the EDM 
constraints, it grows to values of ^ 10 — 100. 

In the scenario proposed here, however, neutralinos are 
not stable, but decay to gravitinos. The resulting grav- 
itino relic density is given in Fig. [21 In the dark green 
shaded region, Q,Qh? is in the range required to account 
for all of non-baryonic dark matter. We see that parts of 
this shaded region are consistent with low energy flavor 
and CP constraints. In this scenario, very large neu- 
tralino freezeout densities are a virtue, not a problem, as 
they allow light gravitinos to have the required relic den- 



FIG. 2: Contours of Q,Qh^ in the (mg,,A) plane,. The thick 
contour is the 2a allowed region. Low-energy constraints and 
fixed GMSB parameters are as in Fig. [1] 
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FIG. 3: Contours of Afs (solid) and lifetime r(x — > G) (dot- 
ted) in the (mg,. A) plane, for A'^m ~ 1, tan/3 = 10, fi> 0, and 
top quark mass mt = 175 GeV. In the light yellow (medium 
blue) shaded region, hadronic (electromagnetic) showers from 
X decays produce discrepancies with BBN observations. The 



sity, despite the significant dilution factor m^/m^. In ^^^^ ^.^^ ^^^^^^^ ^^^2 ^ ^-^ ^ ^j^^ neutralino 



this simple example of minimal GMSB, the Goldilocks 
window, in which both relic density and low energy con- 
straints are satisfied, has uiq ~ 1 — 10 GeV. 

So far, we have considered constraints from parti- 
cle physics and Oqm- We now turn to astrophysi- 
cal constraints. In the preferred band, the gravitino 
is light and dominantly couples through its Goldstino 
components. The neutralino decay widths are T(x ~* 
7G) = {cos^ Bw/^Stt) {my m^M^) and T{x ^ ZG) = 

{sin^ ew/i8TT){ml/m'^^M^) [l - (ml/m^)]^. As shown 
in Fig. [HI these imply lifetimes t > 0.01 s in the pre- 
ferred band. Such late decays are constrained by entropy 
production, /i distortions of the cosmic microwave back- 
ground. Big Bang nucleosynthesis (BBN) fl^ , and small 



LSP region and fixed GMSB parameters are as in Fig. [T] 



scale structure [20]. We find that the last two are most 
stringent, and so focus on them here. 

Standard BBN agrees reasonably well with observa- 
tions. This agreement constrains electromagnetic (EM) 
and hadronic energy release in late decays, which may be 
parameterized by = eiBiY^, where i = EM, had, ei is 
the EM/hadronic energy released in each neutralino de- 
cay, Bi is the branching fraction into EM/hadronic com- 
ponents, and = n^/n^'^, where n^'^ = 2({3)T^/t:^. 
We have determined the following the prescription of 
Refs. 2l| and compared them to the constraints given in 



4 



Ref. [24I. The BBN constraints are shown in Fig. [3] and 
are stringent — in this scenario, neutrahnos are long- 
hved and greatly overproduced, resulting in large energy 
release. In the region of parameter space with 0.097 < 
ilfjh^ < 0.113, the EM (hadronic) constraint requires 
lifetimes r < 10^ s (0.1 s) and > 200 GeV (1 TeV). 

Dark matter produced in late decays also may suppress 
structure on small scales [13]. The free-streaming scale 
Aps — J^^^[''j{t)/a{t)]dt is well approximated by 



Afs - 1-0 Mpc 



r 2 1 


2 


106 S 





1 - 0.07 In 



/ 2 



VlQS s 



(9) 



in the present context, where Ur = IpqI/uIq at decay 
time T, and we have neglected the effect of mz on kine- 
matics and other small effects. Values of Aps are given in 
Fig. [21 they are essentially independent of ttIq. Current 
constraints ^ require Aps ^ 0.2 Mpc, but values near 
this bound may be preferred by observations. Remark- 
ably, constraints from small scale structure are satisfied 
in the region of parameter space allowed by BBN, flavor 
and CP bounds, but just barely — Goldilocks supersym- 
metry therefore predicts "warm" dark matter and may 
explain the suppression of power on scales ~ 0.1 Mpc. 

In summary, we have proposed a simple model in which 
the flavor and overdensity problems of neutralino dark 
matter are simultaneously solved. In the specific frame- 
work considered here, the preferred model is high-scale 
GMSB, with ~ 1 GeV, Vf ~ 10^ GeV, Uy, ~ 100, 
and ~ 2 TeV. This last mass scale is unnaturally 
high, but is dictated by EDM constraints, irrespective of 
cosmology. More generally, this scenario de-emphasizes 
"cosmologically preferred" models with fl^ ~ 0.1, and 
implies that models typically considered excluded by neu- 
tralino overclosure may, in fact, be viable and preferred. 
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